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Abstract 

At strong magnetic fields double-layer two-dimensional-electron-gas systems 
can form an unusual broken symmetry state with spontaneous inter-layer 
phase coherence. The system can be mapped to an equivalent system of 
pseudospin 1/2 particles with pseudospin-dependent interactions and easy- 
plane magnetic order. In this paper we discuss how the presence of a weak 
interlayer tunneling term alters the properties of double-layer systems when 
the broken symmetry is present. We use the energy functional and equations 
of motion derived earlier to evaluate the zero-temperature response functions 
of the double-layer system and use our results to discuss analogies between 
this system and Josephson-coupled superconducting films. We also present 
a qualitative picture of the low-energy charged excitations of this system. 
We show that parallel fields induce a highly collective phase transition to an 
incommensurate state with broken translational symmetry. 

75.10.-b, 73.20.Dx, 64.60.Cn 
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I. INTRODUCTION 



The study of correlated electron systems in fewer than three dimensions continues to 
be an important theme in condensed matter physics. In particular, the study of strongly 
correlated two-dimensional (2D) electron systems on a lattice has been motivated by high- 
temperature superconductivity, while the study of strong correlations in continuum two- 
dimensional systems has been motivated by the fractional quantum Hall Properties 
of high-temperature superconductors are thought by some to be strongly influenced by the 
weak coupling which exists between the superconducting planes. In the fractional quantum 
Hall effect, early work by Halperini anticipated novel fractional quantum Hall effects due 
to interlayer correlationsu in multi-layer systems. Recent technological progress has made 
it possible to produce double-layer two dimensional electron gas systems of extremely high 
mobility in which these effects can be observed. The two electron layers can either be bound 
in separate quantum wells! as illustrated schematically in Fig. (|T|) or bound to opposite 
edges of a single wide quantum well.! In both cases the 2D electron gases are separated 
by a distance d small enough {d ~ 100A) to be comparable to the typical spacing between 
electrons in the same layer. 

In a large magnetic field, strong Coulomb correlations between the layers have long! 
been expected to lead to novel fractional quantum Hall effects. Correlations are especially 
important in the strong magnetic field regime because all electrons can be accommodated 
within the lowest Landau level and execute cyclotron orbits with a degenerate kinetic energy. 
The fractional quantum Hall effect occurs when the system has a gap for making charged 
excitations, i.e. when the system is incompressible. Theory has predicted!!!! that at some 
Landau level filling factors, gaps occur in double-layer systems only if interlayer interactions 
are sufficiently strong. These theoretical predictions have been confirmed experimentally!. 
More recently, theoretical work from several different points of vie has suggested that 
inter-layer correlations can also lead to unusual broken symmetry states with spontaneous 
phase coherence between layers which are isolated from each other (except for inter-layer 
Coulomb interactions). We have arguedS that it is spontaneous interlayer phase coherence 
which is responsible for the recently discovered! extreme sensitivity of the fractional quantum 
Hall effect at total Landau level filling factor v = 1 to small tilts of the magnetic field away 
from the normal to the layers, (y = N/N^ where N is the number of electrons and 
is the number of single-particle levels per Landau level.) In a previous lengthy papeio 
(hereafter referred to as I) we have developed a rather complete description of the novel 
physics associated with spontaneous inter-layer phase coherence in the case where there is 
no tunneling between the layers. The present companion paper will analyze the case of finite 
tunneling between the layers and the role of magnetic field tilt. Many of the ideas presented 



in detail here are discussed qualitatively in Ref. [ ITT 



II. EXPERIMENTAL BACKGROUND 

In this section we review the experimentallll indications that the system is spontaneously 
ordered and exhibits excitations which are highly collective in nature. We focus here and 
throughout this paper on the case of Landau level filling factor v = 1 (that is, 1/2 in each 
layer) and assume for simplicity that the electronic spin degrees of freedom are frozen out 
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by the Zeeman energy. The schematic energy level diagram for the growth-direction degree- 
of-freedom in the double-layer system is shown in Fig. (P for the case of non-interacting 
electrons. For simplicity we assume that electrons can occupy only the lowest electric sub- 
band in each quantum well. If the barrier between the wells is not too strong, tunneling 
from one side to the other is allowed. The lowest energy eigenstates split into symmetric and 
antisymmetric combinations separated by an energy gap Asas which can, depending on the 
sample, vary from essentially zero to many hundreds of Kelvins. The splitting can therefore 
be much less than or greater than the interlayer interaction energy scale, E c = e 2 /ed. 

When the layers are widely separated, there will be no correlations between them and 
we expect no dissipationless quantum Hall state, since each layer hasEl v — 1/2. For smaller 
separations, it is observed experimentally that there is an excitation gap and a quantized 
Hall plateau.Bii This has either a trivial or a highly non-trivial explanation, depending on 
the ratio Asas/E c . For large Asas the electrons tunnel back and forth so rapidly that it is 
as if there is only a single quantum well. The tunnel splitting A S as is then analogous to the 
electric subband splitting in a (wide) single well. All symmetric states are occupied and all 
antisymmetric states are empty and we simply have the ordinary v — 1 integer Hall effect. 
Correlations are irrelevant in this limit and the excitation gap is close to the single-particle 
gap Asas ( or whichever is smaller). What is highly non-trivial about this system is the 
fact that the v = 1 quantum Hall plateau survives even when the tunnel splitting becomes 
arbitrarily small: Asas ^ E c . In this limit the excitation gap has clearly changed to become 
highly collective in nature since the observeoSBi gap can be on the scale of 20K even when 
Asas ~ IK. As we will see below, because of a spontaneous broken symmetry, the 
excitation gap actually survives the limit Asas — ► as illustrated in Fig. (fj). This cross- 
over from single-particle to collective gap is, as we will show, quite analogous to the result 
that for spin polarized single layers, the excitation gap survives the limit of zero Lande g 
factor and hence 'z/ = 1 is a fraction too.'S 

A second indication of the highly collective nature of the excitations can be seen in the 
Arrhenius plots showing thermally activated dissipationE The low temperature activation 
energy A is, as already noted, much larger than Asas- If A were nevertheless somehow 
a single-particle gap, one would expect the Arrhenius law to be valid up to temperatures 
of order A. Instead one observes a rather abrupt leveling off in the dissipation as the 
temperature increases past values as low as ~ 0.1 A. This effect is observed both in double- 
well systems and wide single-well systems and is consistent with the notion of a thermally 
induced collapse of the order that had been producing the collective gap. 

The third significant feature of the experimental data pointing to a highly-ordered collec- 
tive state is the strong response of the system to relatively weak magnetic fields B\\ applied 
in the plane of the 2D electron gases. Within a model that neglects higher electric subbands, 
we can treat the electron gases as strictly two-dimensional. (There is ample evidenced that 
parallel field effects due to subband mixing within a single quantum well produce only small, 
albeit measurable, effects.) B» can then affect the system only if there are processes involv- 
ing tunneling that carry electrons around closed loops containing flux. A prototypical such 
process is illustrated in Fig. (Q). An electron tunnels from one layer to the other at point A, 
and travels to point B. Then it (or another indistinguishable electron) tunnels back and re- 
turns to the starting point. The parallel field contributes to the quantum amplitude for this 
process (in the 2D gas limit) a gauge-invariant Aharonov-Bohm phase factor exp (27r2$/$o) 
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where $ is the enclosed flux and $o is the quantum of flux. Such loop paths evidently con- 
tribute significantly to correlations in the system since the activation energy gap is observed 
to decrease very rapidly with B», falling by factors of order two to ten (depending on the 
sample) until a critical field, £?jj ~ 0.8T, is reached at which point the gap essentially ceases 
changing. To understand how remarkably small B» is, consider the following. We can define 
a length £y from the size of the loop needed to enclose one quantum of flux: £\\Bud = $o- 
(f||[A] = 4.137 x 10 5 /d[A]Sj|'[T].) For Sjf = 0.8T and d = 150 A, £\\ = 0.27//m which is 
approximately twenty times the spacing between electrons in a given layer and thirty times 
larger than the quantized cyclotron orbit radius £ = (hc/eB^) 1 / 2 within an individual layer. 
Significant drops in the excitation gap are already seen at fields of 0.1T implying enormous 
phase coherent correlation lengths must exist. Again this shows the highly collective nature 
of the ordering in this system associated with spontaneous interlayer phase coherence. 



III. SPONTANEOUS PHASE COHERENCE AND SYMMETRY-BREAKING BY 

TUNNELING 

The essential physics of spontaneous inter-layer phase coherence can be addressed either 
from a microscopic point of vie or from a macroscopic Chern-Simons field theory point 

of viewjli'0 but it is perhaps most easily explained in terms of the simple variational wave 
function 

i^)=n{4 T +^4 i }io), (3.i) 

x 

where X is a state label (for instance, the Landau gauge orbital guiding center0) and we are 
using a pseudospin notation in which spin up refers to electrons in the upper layer and spin 
down refers to electrons in the lower layer .0 The interpretation of this wave function is that 
every Landau orbital X is occupied (hence v = 1), but the system is in a coherent linear 
combination of pseudospin up and down states determined by the phase angle (p. This means 
that the system has a definite total number of particles (y — 1 exactly) but an indefinite 
number of particles in each layer. In the absence of inter-layer tunneling, the particle number 
in each layer is a good quantum number. Hence this state has a spontaneously broken 
symmetry0'ii0 in the same sense that the BCS state for a superconductor has indefinite 
(total) particle number but a definite phase relationship between states of different particle 
number. 

In the absence of tunneling, the energy can not depend on the phase angle p and the 
system exhibits a global U(l) symmetry associated with conservation of particle number 
in each layer. One can imagine allowing p> to vary slowly with position to produce excited 
states. Because of the U(l) symmetry, the effective Hartree-Fock energy functional for these 
states has a gradient expansion whose leading term must have the form 

H = ^p s Jd 2 r\Vp\ 2 + ... . (3.2) 

The origin of the finite 'spin stiffness' p s is the loss of exchange energy which occurs when 
ip varies with position. Imagine that two particles approach each other. They are in a 
linear superposition of states in each of the layers (even though there is no tunneling!). If 
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they are characterized by the same phase tp, then the wave function is symmetric under 
pseudospin exchange and so the spatial wave function is antisymmetric and must vanish 
as the particles approach each other. States with spontaneous phase coherence have better 
interlayer electronic correlations and hence lower interlayer Coulomb interaction energy. If 
a phase gradient exists then there is a larger amplitude for the particles in opposite layers 
to be near each other and hence the interlayer interaction energy is higher.E3 This loss of 
exchange energy is the source of the finite spin stiffness and is what causes the system to 
spontaneously 'magnetize'. 

A skeptical reader might legitimately worry that the single-Slater-determinant varia- 
tional wavefunctions in terms of which we have framed the above discussion misrepresent 
the physics. Indeed, the ground state of double-layer systems in the Hartree-Fock approx- 
imation has spontaneous phase coherence at all values of Vt whereas we believe that this 
property actually holds only at a discrete set of total filling factors, including v T — 1. To 
understand why these concerns would be misplaced, it is useful to briefly review some of 
the microscopic physics underlying the incompressible states whose occurrence is responsi- 
ble for the fractional quantum Hall effect. For a pair of interacting electrons confined to 
the lowest Landau level, only one relative motion state is available for each relative angular 
momentum. This property means that the interaction is characterized by a discrete set of 
energy scales, VJ, first identified by Haldane@ and known as Haldane pseudopotentials. V\ 
is the interaction energy of a pair of electrons with relative angular momentum /; in the 
double-layer case, the pseudopotential values for interlayer and intralayer interactions will 
differ. Many of the largest charge gaps which occur in the fractional quantum Hall regime 
can be understood in terms of Haldane pseudopotentials. For example, for double-layer 
systems with nearby layers the largest energy scales should be the interlayer and intralayer 
/ = pseudopotentials. It is easy to show that it is possible to form many body-states in 
which pairs of particles never have I = for < 1. This is not possible for > 1. Thus 
if Vf — 1, the energy to add a charge to the system is greater than the energy to remove 
one, hence the charge gap. For — 1 it is easyil to show that the only states which com- 
pletely avoid pairs with relative angular momentum I = are, up to rotation in pseudospin 
space, identical to the wavefunctions in Eq. ( |3.1| ). This is the real reason why the variational 
wavefunctions used above are accurate. In the following we describe low-energy states of the 
system as spin-texture states with a position-dependent local pseudospin orientation. The 
property that the low-energy states are completely specified by spin-textures, which we will 
use frequently in following sections, depends on both the broken symmetry of the ground 
state and on the existence of the charge gap. 

The U(l) symmetry leads to Eq. fl3.2| ) which defines an effective XY model which will con- 
tain vortex excitations which interact logarithmically. In a thin film of superfluid 4 He, vor- 
tices interact logarithmically because of the energy cost of supercurrents circulating around 
the the vortex centers. (In superconducting thin films the same logarithmic interaction 
appears but is cut off on length scales exceeding the penetration depth.) Here the same log- 
arithmic interaction appears. Microscopically this interaction is due to the potential energy 
cost (loss of exchange) associated with the phase gradients (circulating pseudospin currents). 
Hartree-Fock estimates^ indicate that p s and hence the Kosterlitz-Thouless critical temper- 
ature are on the scale of IK in typical samples. Vortices in the tp field ('merons'B) are 
reminiscent of Laughlin's fractionally charged quasiparticles but in this case carry charges 



5 



±~e and can independently be left- or right-handed for a total of four 'flavors'. Bound 
meron pairs with opposite vorticity are the lowesti! energy charged excitations of the sys- 
tem. The finite pseudospin stiffness not only permits the presence of spontaneous pseudospin 
magnetization but leads to a finite charge excitation gap (even though the tunnel splitting 
is zero). Thus the QHE survivesB@ the limit Asas — ► 0. 

Since the 'charge' conjugate to the phase <p is the z component of the pseudo spin S z , 
the pseudospin 'super cur rent' 



J 



2p* 

h 



(3.3) 



represents oppositely directed charge currents in each layer. Below the KT transition tem- 
perature, such current flow will be dissipationless (in linear response) just as in an ordinary 
superfluid. Likewise there will be a linearly dispersing collective Goldstone mode as in a 
superfluid.SSSIffl 

To reinforce the idea that this is not an ordinary superfluid or superconductor, it is 
perhaps useful to rewrite the original variational wave function as 



|V>T>, 



where 



(3.4) 



(3.5) 



x 



is the fully up polarized spin state. We now see that the analogy to an excitonic insulatoiili! 
is closer than the analogy to a superconductor. Only a gauge-neutral object (a particle bound 
to a hole) can condense and propagate freely in a strong B field. A phase gradient Vy? causes 
a flow of these neutral objects which then constitutes a spin current analogous to the charge 
current in a superconductor. In analogy to the excitonic insulator case, one can make a 
particle-hole transformation in one of the layers to produce a formal resemblance to the 
BCS mean field state: 



n 



X I 



1 

V2 



(1 + 



e%<P c x,i c x,\ 



l^t>, 



(3.6) 



but it is then clear that one has a pairing between a particle and holejnot between two parti- 
cles. It is worth remarking that there is in general an exact mappingEa between double-layer 
electron-hole systems and two-component electron systems in the limit of strong magnetic 
fields. The Hamiltonian for a double-layer electron system is mapped, up to a constant 
term, to that for an electron-hole system if a particle-hole transformation is made in one of 
the layers. The filling factor of interest here, vt — 1, corresponds to an electron-hole sys- 
tem with equal electron and hole densities and, as explained above, the spontaneous phase 
coherence state corresponds to an excitonic superfluid state. In the electron-hole case, the 
fact that a broken symmetry ground state can occur had been appreciated some years ago.il 
Some of the recent advances in understanding the physics associated with the spontaneous- 
phase-coherence ground state in the electron-electron case, have important implications for 
the electron-hole case which we will not explore at length here. 
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A finite tunneling amplitude t between the layers breaks the U(l) symmetry 



#nff = d\ 



(3.7) 



by giving a preference to symmetric tunneling states. This can be seen from the tunneling 
Hamiltonian 



H T = -tj d 2 r {^(r)Vi(r) +^(r)^ T (r)} (3.8) 
which can be written in the spin representation as 

H T = -2t f d 2 rS x (r). (3.9) 



(Recall that the eigenstates of S x are symmetric and antisymmetric combinations of up and 
down.) 

We can shed further light on the spontaneous symmetry breaking by considering the 



tunneling Hamiltonian in Eq. (3.8) as a weak perturbation. Naively, since particle 



number is separately conserved in each layer for t = 0, one might expect 

lhn^(iP\H T \if;) = 0. (3.10) 

That is, one might expect that the first-order term in the perturbation series for the energy 
due to t to vanish. Instead however we find that the energy shifts linearly in t 



lim lim — (if) \H^\ if>) = lim lim —rii^ 



d 2 r2SJv 



■0 



t ►OA >oo tA t — >0i >oc A 

= -m x , (3.11) 

where A is the system area, and m x is, by definition, the magnetization which is the system's 
order parameter.!! If the interlayer spacing d is taken to be zero, one can readily showlil 

that the variational wave function in Eq. (|3.1| ) is exact, hence linn, ^m x = 1, and t = 

Asas/2- For finite d, Eq. ( |3.1|) is no longer exact and quantum fluctuations will0 reduce 
the magnitude of m x and we must renormalize the hopping parameter t appropriately. 

As the layer separation d increases, a critical point d* will be reached at which the 
magnetization vanishes and the ordered phase is destroyed by quantum fluctuations. 00 This 
is illustrated in Fig. (Q). For finite tunneling t, we will see below that the collective mode 
becomes massive and quantum fluctuations will be less severe. Hence the phase boundary 
in Fig. (0) curves upward with increasing A S as- For A S as = the destruction of long-range 
order and the charge excitation gap are intimately related and occur simultaneously at d* 
and zero temperature. For finite As as the system always has non-zero m x even in the phase 
with zero charge gap. 

The effective Hamiltonian in Eq. ( |3.7| ) looks like the sine-Gordon model which is known 
to have a finite temperature Kosterlitz-Thouless phase transition.i3 One might be tempted 
to speculate then that the rapid collapse of the Arrhenius plots of the dissipation at unex- 
pectedly low temperatures is associated with a true phase transition. We believe however 
that these are merely rapid cross-overs rather than true phase transitions because the phase 
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ip is compact. The quantum states denned by ip and tp + 2ir are microscopically identical. 
The interpretation of ip in a sine-Gordon theory as an electrostatic potential for Coulomb 
charges or as the surface height in the solid-on-solid model requires that ip and ip + 2tt be 
distinguishable states. Hence what we really have is an XY model in a symmetry breaking 
field which has no true phase transition since the vortices are linearly confined. [Note that 
the solid-on-solid model has no analog of open strings (see Sec. |TV| ) terminated by vortices 
such as we have here.E-J| 



IV. CHARGED EXCITATIONS 



At filling factor v — 1, there is an intimate connection between local distortions of the 
pseudospin orientation and the local charge density.0 A simple example of this is provided 
by the fully spin polarized v = 1 state of Eq. ( |3.1|) . Since the Landau level is filled, the 
charge density is uniform and the Pauli principle forbids any (intra-Landau-level) excitations 
which do not flip spins. One can form a spin-flip particle-hole pair at an energy cost o 



2 , , I 

e /7r\ 2 



u > = AsAS + Tl UJ (4 - 1) 

where the first term is the tunneling energy necessary to flip the pseudospin and the second 
represents the loss of Coulomb exchange energy (here computed for the special case d = 0). 
The Coulomb exchange energy gives a finite cost even in the absence of tunneling and 
this explains the fact that the typical observed gap can be much larger than the tunnel 
splitting]!!'! The simple spin flip is not however the optimal charged-pair excitation. From 
analogy with results of Sondhi et al.0 For the case of 'real' spins, we know that smooth local 
distortions of the pseudospins produce a charge density given by the remarkable formula§0 

6p(r) = -—e iat m(r) • [<9 M m(r) x d u m(r)] (4.2) 

where m is a unit vector giving the local pseudospin orientation. Sp(r) is exactly v times the 
Pontryagin index, or topological charge density.0'0 The density in Eq. (|4.2| ) can be viewed 
as the time-like component of a conserved (divergenceless) topological 'three-current' 

J a = ~ e a ^ e abc m a {v)d p m b {v)d 1 m c {v). (4.3) 
on 

Using the fact that m is a unit vector, it is straightforward to verify that <9^j M = 0. 

'Skyrmion' (hedgehog) configurations of the order parameter carry net charge ±1 and, 
when Asas = and d = 0,00 a skyrmion pair has an excitation energy which is one-half 
of that of a simple spin-flip pair for the case of Coulomb interactions between the electrons. 
At finite d, the SU(2) symmetry is lowered to U(l) and the cheapest charged excitations are 
composed of 'merons' which are essentially vortex solutions in which the local pseudospin 
winds by ±2n at infinity and tilts either up or down out of the XY plane in the core region 
as shown in Fig. (9) of I. Integration of the charge density using Eq. (|4.2| ) shows that vortices 
carry charge ±|. They are somewhat analogous to Laughlin quasiparticles, however they 
differ considerably in that, below the Kosterlitz-Thouless temperature, they are confined 
together in vorticity neutral pairs by their logarithmic interaction. The cheapest object 



S 



with a net charge is then a vortex-antivortex pair, with each vortex carrying charge +| 
(or — ~) for a total charge of +1 (or —1). The charge excitation cost can be estimated by 
minimizing 



E- 



pair 



2E mc ~F' 



AeR 



2iip s In 



R 



(4.4) 



where E mc is the meron core energyil , and -R mc is the meron core size. The optimal 
separation is given byE-a Rq = e 2 / (8nep s ). It is important for the discussion below that in 
typical double-layer systems p s is much smaller than the microscopic energy scale e 2 'jet. 
For djl — 1, p s — 6.19 x 10 _3 (e 2 jet) in the Hartree-Fock approximation and it is further 
renormalized downward by quantum fluctuations.^ Typical values of p s for double-layer 
systems are smaller than 5 x 10~ 3 (e 2 jet) so that Ro/t will typically be larger than ~ 8t. 
The small values of the pseudospin stiffness allow the charged pseudospin textures to be 
large, as required for the validity of the long-wavelength description being employed here. 

The introduction of finite tunneling amplitude destroys the U(l) symmetry and makes 
the simple vortex-pair configuration extremely expensive. To lower the energy the system 
distorts the spin deviations into a domain wall or 'string' connecting the vortex cores as 
shown in Fig. (f|). The spins are oriented in the x direction everywhere except in the domain 
line region where they tumble rapidly through 2ir. The domain line has a fixed energy per 
unit length and so the vortices are now confined by a linear potential corresponding to a 
fixed 'string tension' rather than being confined only logarithmically. We can estimate the 
string tension by examining the energy of a domain line of infinite length. The optimal form 
for a domain line lying along the y axis is given by 



(f(r) = 2arcsin [tanh 
where the characteristic width of the string is 

'27Tt 2 Ps 



(4.5) 



(4.6) 



The resulting string tension ist 



T = 8 



tps 

2ixi 



5ps 



(4.7) 



Provided the string is long enough (R ^> £), the total energy of a segment of length R will 
be well-approximated by the expression 



EL 



pair 



2E'. + 



AeR 



+ T R. 



(4.8) 



The prime on E mc in Eq. ( f4.8|) indicates that the meron core energy can depend on As^s- 
is minimized at R = R' Q = we 2 /4eXo where it has the value 



EL 



pair 



E 



pair 



2E' mc + Je 2 T /e. 



(4.9) 
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Note that apart from the core energies, the charge gap at fixed layer separation (and hence 
fixed p s ) is oc T^ 2 oc t 1 / 4 ~ ^sas^ which contrasts with the case of free electrons, for which 
the charge gap is linearly proportional to Asas- Note that because the exponent 1/4 is so 
small, there is an extremely rapid initial increase in the charge gap as tunneling is turned 
on. 

The crossover between the meron-pair pseudospin texture which holds for t = and the 
domain line string pseudospin texture described above occurs at a finite value of t which 
we can estimate by the following argument. For R' > Rq the vortices are already bound 
by the logarithmic attraction due to the gradient energy before the linear attraction due 
to the hopping becomes important at larger separations. In this regime tunneling does not 
play an important role in determining the nature of the lowest energy charged pseudospin 
texture. As t increases R' oc t -1 / 4 decreases and will eventually reach Rq which is, of course, 
independent of t. Since 

the characteristic width of the domain line becomes comparable to R' in the same range 
of t values where R' and Rq become comparable. We may conclude that the nature of 
the charged pseudospin texture crosses over directly from the meron pair form to the finite 
length domain line string form for p s /(e 2 /e£) ~ 1/25, or equivalently for t ~ t CT where 



4 x 10 d 



e 2 /e£ 



3 e 2 



(4.1D 



The crossover tunneling amplitude is thus typically smaller than 5 x 10~ 4 (e 2 /e£). Typical 
tunneling amplitudes in double-layer systems are smaller than ~ 10 -1 (e 2 /e£) and can be 
made quite small by adjusting the barrier material or making the barrier wider. Never- 
theless, it seems likely that t will be larger than t cr except for samples which are carefully 
prepared to make t as small as possible. As t increases beyond t CT , R' will continue to 
decrease. When R' becomes comparable to the microscopic length £, the description given 
here will become invalid and the lowest energy charged excitations will have single-particle 
character. However, the domain-wall string picture of the charged pseudospin texture has a 
very large range of validity since R' oc t -1 / 4 decreases very slowly with increasing t. Writing 
R ~ (e 2 /8enp s )(t CI /t) 1 / 4 we find that R' ~ I only for t ~ l(T 2 [(e 2 /el) 2 / p s ). Using typical 
values of p s we see that the charged excitation crosses over to single particle character only 
when the hopping energy t becomes comparable to the microscopic interaction energy scale. 
The various regimes for the charge excitations of double-layer systems are summarized in Ta- 
ble |. Almost all typical double-layer systems lie within the regime of the domain- wall-string 
pseudospin texture charge excitation. 

We should emphasize that all the discussion of charged excitations above assumes that 
the meron core energies do not make a dominant contribution to the charge excitation en- 
ergies and that meron core sizes are small compared to the overall size of the quasiparticles. 
Recent calculations!! have demonstrated that these conditions are never well satisfied when 
the Hartree-Fock approximation is used to approximate the charged excitations. The phys- 
ical pictures summarized still have some qualitative validity, however. The Hartree-Fock 
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approximation neglects quantum fluctuation effects which reduce both the spin- stiffness and 
the order parameter, tending in both cases to increase the size of the quasiparticles and 
increase the appropriateness of the pictures presented here. 



Murphy et al.l-tta and Santos et alB have shown that the charge gap in double layer 
systems is remarkably sensitive to the application of relatively weak magnetic fields B», 
oriented in the plane of the 2D electron gas. Experimentally this field component is generated 
by slightly tilting the sample relative to the magnetic field orientation. Tilting the field (or 
sample) has traditionally been an effective method for identifying effects due to (real) spins 
because orbital motion in a single-layer 2DEG system is primarilyc3 sensitive to B±, while 
the (real) spin Zeeman splitting is proportional to the full magnitude of B. Adding a parallel 
field component will tend to favor more strongly spin-polarized states. For the case of the 
double layer v = 1 systems studied by Murphy et al.Jl the ground state is known to already 
be an isotropic ferromagnetic state of the true spins and the addition of a parallel field would 
not, at first glance, be expected to influence the low energy states since they are already 
fully spin-polarized. (At a fixed Landau level filling factor B± is fixed and so both the total 
B and the corresponding Zeeman energy increase with tilt). Nevertheless experiments! have 
shown that these systems are very sensitive to B\\ . The activation energy drops rapidly (by 
factors varying from two up to an order-of-magnitude in different samples) with increasing 
B\\. At B\\ = Bt there appears to be a phase transition to a new state whose activation gap 
is approximately independent of further increases in B». 

The effect of B\\ on the pseudospin system can be visualized in two different pictures. In 
the first picture we use a gauge in which By = V x Ay where An = B\\ (0, 0, x). In this gauge 
the vector potential points in the z direction (perpendicular to the layers) and varies with 
position x as one moves parallel to the layers. As an electron tunnels from one layer to the 
other it moves along the direction in which the vector potential points and so the tunneling 
matrix element acquires a position-dependent phase t —> t e l ® x where Q = 2n / 'Ln and 
L\\ = $> /Bnd is the length associated with one flux quantum $ between the layers [defined 
in Fig. (^|)]. This modifies the tunneling Hamiltonian to Ht = — J d 2 r h(r) ■ S(r) where h(r) 
'tumbles': i.e., h(r) = 2t (cosQx, sin Qx, 0). The effective XY model now becomes 



which is precisely the Pokrovsky-Talapov (P-T) modelEj and has a very rich phase diagram. 
For small Q and/or small p s the phase obeys (at low temperatures) ip(r) = Qx; the moment 
rotates commensurately with the pseudospin Zeeman field. However, as Bn is increased, the 
local field tumbles too rapidly and a continuous phase transition to an incommensurate state 
with broken translation symmetry occurs. This is because at large B\\ it costs too much 
exchange energy to remain commensurate and the system rapidly gives up the tunneling 
energy in order to return to a uniform state Vy? ~ which becomes independent of B\\. 
As explained in further detail below weEl find that the phase transition occurs at zero 
temperature for 
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(5.1) 
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B\ = B± (2£/7rd)(2t/7T Ps ) 1 / 2 . (5.2) 

Using the parameters of the samples of Murphy et al.S and neglecting quantum fluctuation 
renormalizations of both t and p s we find that the critical field for the transition is ~ 1.6T 
which is within a factor of two of the observed valuefi Note that the observed value Bu = 
0.8T corresponds in these samples to a large value for L»: L»/£ ~ 20 indicating that the 
transition is highly collective in nature. We emphasize again that these very large length 
scales are possible in a magnetic field only because of the interlayer phase coherence in the 
system associated with condensation of a neutral object. 

Having argued for the existence of the commensurate-incommensurate transition, we 
must now connect it to the experimentally observed transport properties. In the commen- 
surate phase, the order parameter tumbles more and more rapidly as B\\ increases. As we 
shall see below, it is this tumbling which causes the charge gap to drop rapidly. In the 
incommensurate phase the state of the system is approximately independent of B\\ and this 
causes the charge excitation gap to saturate at a fixed value. 

Recall that in the presence of tunneling, the cheapest charged excitation was found to be a 
pair of vortices of opposite vorticity and like charge (each having charge ±1/2) connected by 
a domain line with a constant string tension. In the absence of B\\ the energy is independent 
of the orientation of the string. The effect of B\\ is most easily studied by changing variables 
to 

0(r) = <p(r) - Qx. (5.3) 

This variable is a constant in the commensurate phase but not in the incommensurate phase. 
In terms of this new variable, the P-T model energy is 

H = J >r |i p s [(d x 9 + Qf + {d y 9f} - JL cos* j. (5.4) 

We see that B\\ defines a preferred direction in the problem. Domain walls will want to line 
up in the y direction and contain a phase slip of a preferred sign (— 2n for Q > 0) in terms 
of the field 0. Since the extra term induced by Q represents a total derivative, the optimal 
form of the soliton solution is unchanged. However the energy per unit length of the soliton, 
which is the domain line string tension, decreases linearly with Q and hence B»: 

(5.5) 

where T is the tension in the absence of parallel B field given by Eq. ( f4.7| ) and B?, is the 
critical parallel field at which the string tension goes to zero.il We thus see that by tuning 
By one can conveniently control the 'chemical potential' of the domain lines. The domain 
lines condense and the phase transition occurs when the string tension becomes negative. 

Recall that the charge excitation gap is given by the energy of a vortex pair separated 
by the optimal distance Rq = ^/e 2 /4eT. From Eq. ( |4.8|) we have that the energy gap for 
the commensurate state of the phase transition is given by 

i 

e T/e 



T = Tn 
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B \\ 
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(5.6) 

As B\\ increases, the reduced string tension allows the Coulomb repulsion of the two vortices 
to stretch the string and lower the energy. Far on the incommensurate side of the phase 
transition the possibility of interlayer tunneling becomes irrelevant. From the discussion of 
the previous section it follows that the ratio of the charge gap at B\\ = to the charge gap 
at B\\ — > oo should be given approximately by 

^ = (t/t cr ) 1/4 « (e 2 /^) 1/2 t 1/4 8p, 3/4 . (5.7) 

Putting in typical values of t and p s gives gap ratios ~ 1.5 — 7 in agreement with experiment. 
According to the discussion of the previous section, gap ratios as large 
0.07(e 2 /e£)/p s , can be expected in the regime where the pseudospin texture picture applies. 
Here t max is the hopping parameter at which the crossover to single-particle excitations 
occurs. Thus gap ratios as large as an order of magnitude are easily possible. Of course, 
all the discussion here neglects orbital effects (electric subband mixing) within each of the 
electron gas layers, and these will always become important at sufficiently strong parallel 
fields. 

It should be emphasized that only this highly collective picture involving large length 
scale distortions of topological defects can possibly explain the extreme sensitivity of the 
charge gap to small tilts of the B field. Recall that at BT, the tumbling length L\\ is much 
larger than the particle spacing and the magnetic length. Simple estimates of the cost to 
make a local one-body type excitation (a pseudospin-flip pair for example) shows that the 
energy decrease due to B\\ is extremely small since £/L\\ is so small. As we will see in Sec. 
VI I| numerical exact-diagonalization calculations on small systems confirm the existence of 
this phase transition and show that the fermionic excitation gap drops to a much smaller 
value in the incommensurate phase. 

We now discuss the commensurate-incommensurate phase transition from the micro- 
scopic point of view. At d — 0, the B\\ = Landau-gauge many-body ground state wave- 
function is a single Slater determinant in which the single-body states are the symmetric 
linear combination of two single-layer states with the same guiding center. This is the state 
represented by Eq. ( |3.1| ) with (p = 0. Phase coherence is established (either spontaneously 
or in this case) by tunneling between single-layer states with the same guiding center. For 
many purposes this state is still a good approximation to the ground state at finite d since 
it optimizes the tunneling energy and has good correlation energy; an electron in one layer 
automatically sees an exchange-correlation hole in the other layer at the same place. (It 
would remain an exact ground state at finite d in the absence of interactions.) From a mi- 
croscopic point of view it is the good interlayer correlations of states with phase coherence 
which leads to the broken symmetry in the absence of tunneling. 

To see the effect of a parallel B field, it is convenient to choose a new Landau gauge 
for the perpendicular field = — B±(y, 0, 0). In this gauge, a parallel field giving rise to 
t — > texp(iQx) causes tunneling to couple states in the two layers whose wavevectors differ 
by Q and whose guiding centers therefore differ by Q^ 2 .H Thus, for non-interacting electrons 
the exact ground state in a parallel field is one in which the exchange-correlation hole is not 



A + Je*T Q /e 
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directly opposite its electron but rather is shifted away by £ 2 Qy as the B field tilts in the x 
direction (i.e., the displacement is perpendicular to the direction of the in-plane field): 

i^^nRt+^ji )' ( 5 - 8 ) 

Y 

This state maintains all of its tunneling energy but rapidly loses interlayer correlation energy 
as the field tilts. At large tilt it is better to give up on the tunneling by shifting the two 
layers relative to each other to put the correlation hole back next to its electron. 

This shift can be seen to be the change from commensurate to incommensurate states 
discussed above. A straightforward computation shows that the commensurate state has 
the pseudospin tumbling 

<^|(r)^(r)|^ Q > = e-i^e-"* (5.9) 

while the pseudospin is constant in the incommensurate phase. 

All of our discussion of the phase transition in a parallel field has been based on mean- 
field theory. Close to the phase transition, thermal fluctuations will be important. At finite 
temperatures there is no strict phase transition at Bt in the the P-T model. However there 
is a finite temperature KT phase transition at a nearby B\\ > Bu. At finite temperatures 
translation symmetry is restoredll in the incommensurate phase by means of dislocations 
in the domain string structure. Thus there are two separate KT transitions in this system, 
one for t — 0, the other for t 7^ and B\\ > Recently Readi has studied this model 
at finite temperatures in some detail and has shown that at the critical value of B» there 
should be a square-root singularity in the charge gap on both sides of the transition. The 
existing data does not have theresolution to show this however. At zero-temperature the 
commensurate- incommensurate^'^ phase transition must be treated quantum mechanically. 
It is necessary to account for the world sheets traced out by the time evolution of the strings 
which fluctuate into existence due to quantum zero-point motion. Read also points out that 
the inevitable random variations in the tunneling amplitude with position, which we have 
not considered at all here, cause a relevant perturbation. 

VI. COLLECTIVE MODES AND RESPONSE FUNCTIONS 

In this section we will discuss the charge neutral collective excitations of double-layer 
systems and some physically important response functions which have poles at the collective 
excitation energies. In the pseudospin language, the collective excitations are spin-waves in 
which the pseudospin precesses around its ground state orientation. We will thus need to 
enlarge our description of the system by allowing the pseudospin texture to have orientations 
out of the x — y plane. This requires that we generalize from the AF-model description 
of the system employed in previous sections, to the more complete anisotropic non-linear 
sigma model description which we have discussed at length in Ref. [ [TJj (I). (Actually this 
generalization is also required if want to render the physics of the meron cores.) States of the 
system are characterized by a pseudospin texture function, m(r) which specifies the space 
dependence of the pseudospin orientation. The energy of a pseudospin texture is given by 
the following functional (where we retain only the leading terms in number of derivatives): 
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E[m(r)] 



d\ 



/3(m 2 z (r)) + -^-|Vm_L| 2 — nt m a 



(6.1) 



where m(r) • m(r) = 1 and m±_ is the projection of the unit vector onto the x — y plane. The 
energy is with respect to the ground state in the absence of tunneling. Here n = (27r£ 2 ) -1 
is the total electron density. This energy functional is expected to be accurate for states 
with spin-textures which vary slowly on a microscopic scale. If the pseudospin is confined to 
the x — y plane, its orientation is specified by its azimuthal angle and the energy functional 
reduces to the XY functional used above. This Ginzburg-Landau energy functional differs 



from the one discussed in Ref . [ [14]] only through the addition of the tunneling term. The form 
of the functional follows from symmetry considerations and the parameters can be considered 

I4}| , Eq. ( |6.1| ) was derived microscopically using a 



as phenomenological constants. In Ref. [ 
Hartree-Fock approximation, and explicit expressions for the coefficients were obtained, 
which become exact in the limit, d — > 0. We now generalize the discussion of collective 
modes given there to the case where tunneling occurs between the layers. 



A. Collective Modes with Tunneling: Bn = 

In the presence of inter-well tunneling, the pseudospin-orientation in the ground state is 
constant and points in the x direction. To calculate response functions we use the equations 
of motion derived in Ref. [ ITJ| . We will assume that for low-lying excitations the pseudospin 



orientation is always close to the x direction. Then the equation of motions of the pseudospin 
texture in Fourier space is: 

dm y (q) _ 4tt£ 2 dE[m] ^ 



dt h dm z {—Q 

dm z (q) _ 4n£ 2 dE[m] 
dt h dm y (—q) 



(6.3) 



In the rest of this section we will use £ as the unit of length. It is possible to linearize 
these equations if the pseudospin orientation is close to the x direction by letting m x = 
1 — {m 2 + m 2 z )/2 + • • •, and dropping terms higher than first order in m y and m z . For the 
double-layer system m y is proportional to the current flowing locally between layers and m z 
is proportional to the difference between the local densities in the two layers. To calculate 
the response functions of interest we add terms to the energy functional corresponding to 
pseudospin Zeeman fields in the y and z directions, h y and h z . Physically h y corresponds 
to a perturbation in which an imaginary term is added to the tunneling amplitude and h z 
corresponds to a bias potential between the two layers. Linearizing and adding the Zeeman 
terms we find that 

' 1 - -{2/3 + tn)m z {q) + h z (6.4) 



dt h 

(tn + p s q 2 )m y (q) - h y . (6.5) 



dm z (q) An ,, 

— -j — = ~f[y tnJr ) ,h >i - m ' - ' '/ ■ 



There is a similarity between these equations of motion and those of the Josephson 
effect, which is connected to the similarities between inter layer phase coherence and su- 
perconductivity mentioned previously.^ The current across a Josephson junction between 
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two superconductors is proportional to sin(0) where is the difference in the phase of the 
order parameter across the junction. The current between layers in a double-layer system 
is similarly proportional to m y , i.e. proportional to sin(0) where <fi in the double-layer case 
specifies orientation of the component of the pseudospin in the x — y plane. [In the linearized 
case we are discussing we can equate and sin(</>).] More physically (ft in the double- layer 
case specifies the difference between the phase-coherence angle and the phase angle of the 
interlayer tunneling amplitude. In the case of the Josephson effect 

~dt = — (6 ' 6) 

where V is the potential drop across the junction. Thus the equations of motion for the phase 
angles in the Josephson effect case and in the present case differ because of the presence 
of the term proportional to (2/3 + tn)m z (q)] a Josephson effect could be achieved if the 
tunneling current were somehow extracted from the double-layer system sufficiently quickly 
to prevent any difference in the charge densities of the two-layers from building up, i.e., if 
m z where identically zero. In the case of double-layer systems, unlike the case of a Josephson 
junction between two superconductors, it seems to be impossible to do this. The results we 
derive below are for an isolated double-layer system and the resonance frequencies we obtain 
are analogous to the Josephson plasma oscillations in an isolated Josephson junction. 

Solving for the pseudospin magnetization induced by Zeeman fields with frequency uj we 
obtain the following result for the pseudospin response tensor: 



m y \ h I —iu 87r/3 + 2t\ ( —4iTh z 



v 

m 



(6.7) 



" E 2 cm -(huj) 2 \-{2t + ^q 2 Ps ) -iuj j \ 47rh y J 
where the collective mode energy is given by 

E cm = [(2t + 8n(3)(2t + 4vrg 2 p s )] 1/2 - (6-8) 

The linearly dispersing collective mode of the t — > case acquires a gap because of the lifting 
of the U(l) symmetry. 

All components of this response tensor have poles at the collective mode energies. The 
response of the charge density difference to a time-dependent inter-layer bias potential is 
given by 

_^rn 1 _ (4nh)2t + A-nq 2 p s 
XZZ ~ K~ E 2 cm - (hu) 2 {bJ) 

Using the continuity equation we can evaluate the corresponding conductivity, the response 
in oppositely directed electric currents to oppositely directed electric fields in the two-layers: 

a zz (q,u) = e 2 uxzz(q,u)/iq 2 (6.10) 

For t — > the real part of the conductivity has a 5-function peak at zero frequency leading 
to spin-channel superfluidity.0 In the presence of interlayer tunneling the ^-function peak is 
shifted to u = 2t/h and the superfluid behavior is lost. 

In the static (u — > 0) long- wavelength (q 0) limit Xzz approaches a constant: 
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(8tt/5 + 2t) ' 

A constant static interlayer bias potential will@ simply tilt the pseudospin orientation 
slightly out of the x — y plane. The effect of interlayer tunneling is to favor smaller tilts. On 
the other hand 

Xw (u, = Q,q = 0) = ± (6.12) 

If the tunneling amplitude goes to zero there is no restoring force for rotations of the pseu- 
dospin in the x — y plane, the addition of an infinitesimal imaginary tunneling amplitude will 
shift the pseudospin orientation from the x direction to the y direction, and Xyy{ u = 0, = 0) 
diverges. 



B. Collective Modes with Tunneling: Commensurate State 

For B« = we have been able to calculate collective mode energies and linear response 
functions by linearizing the non-linear sigma model energy functional around the ground 
state pseudospin orientation. For Bn ^ 0, the ground state pseudospin orientation in the 
commensurate state rotates with the pseudospin Zeeman field and it is necessary to linearize 
the non-linear sigma model energy functional around the rotating pseudospin texture. In 
anticipation of our needs for the case of the incommensurate state we allow an arbitrary 
rate of rotation for the rotating frame pseudospin function: 

fh x = m x cos Px — m y sin Px (6.13) 
fhy = m x smPx + rriyCosPx. (6.14) 

In order to be able to treat the problem analytically, we limit ourselves to the case of a 
position-independent tumbling rate P. The Ginzburg-Landau energy functional, expressed 
in terms of the rotating frame pseudospins, is: 

zji f j2 a 2,Ps|v7~ 12 i l ~ 1 2 

E = Jar pm z + — |Vmj_| H — |mj_| 

—tn[m x cos(P — Q)x + m y sin(P — Q)x] 



(6.15) 



If we allowed only translationally invariant spin-textures in the rotating frame for which the 
lowest energy occurs, this energy functional would be minimized for small parallel fields by 
choosing P = Q to obtain an energy per area of p s Q 2 /2 — tn. This is the commensurate 
state. At large parallel fields the energy functional would be minimized by the state with 
P = 0, which has the same energy as if no tunneling were present. These two states cross 
in energy when Q 2 ^ 2 = 2. However, as we detail in the following section, the ground state 
at high parallel fields can lower its energy further by breaking translational symmetry. The 
commensurate state is the ground state only for Q < Q c where Q 2 i 2 = 16 /n 2 . Q = Q c for 

B \\ = Bi- 
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In order for the linearization of this functional to be valid we must assume we can choose 
P so that m x is close to 1 everywhere, both for the ground state and for the collective 
excitations in which we are interested. Assuming this to be the case we can apply periodic 
boundary conditions in the thermodynamic limit and drop the total derivative term in the 
energy functional. For the commensurate state the linearized energy functional simplifies to 
the following Fourier space version: 

E[m y , mJ=^(^-tn) + £(/3 + f- ^) K(q)| 2 + (f g 2 + f ) \m y \ 2 . (6.16) 

The first term on the right-hand-side of Eq. ( |6.16| ) is the ground state energy of the commen- 
surate state. The fact that it becomes positive at large Q implies that the commensurate 
state eventually becomes unstable. As discussed in previous sections we expect the state to 
become unstable to the introduction of phase-slips or solitons, which involve only the planar 
portion of the pseudospin texture, when p s Q 2 > 16tn/ir 2 . We note from Eq. ( |6.16| ) that if (3 
is small enough the ground state will become unstable with respect to textures with m z ^ 
before the solitons are introduced. The requirement for the applicability of the scenario 
we introduced previously based on the XY language is that (3£ 2 > t (16 — n 2 )/ (47r 3 ). From 
estimates of (3 in Ref. [ [14]], it is clear that for present double-layer samples this condition 
is secure. Nevertheless, the possibility exists that the behavior in parallel fields could be 
quite different from that described here, for double-layer systems with a layer spacing much 
smaller than what is achievable at present. We have not explored this regime in detail and 
will assume in what follows that 2(3 + tn — p s Q 2 is positive in the incommensurate state. 
The effect of B\\ is then to reduce the collective mode energy, just as the charged excitation 
energies are lowered. Solving the equations of motion with the linearized energy functional 
of the commensurate state we find that 



E, 



cm 



(2t + 87r/3 - Anp s Q 2 )(2t + Anp s q 2 )} 1/2 (6.17) 



This result agrees with collective mode energies calculated in the Hartree-Fock 
approximation^ if the Hartree-Fock approximation is used for the pseudospin stiffness. 



C. Collective Modes with Tunneling: Incommensurate State 

For B\\ > -By it becomes energetically favorable to break translational symmetry and 
introduce phase slips, or solitons, into the pseudospin texture. In the ground state the 
phase slips are periodic00 with a period L s which is determined by minimizing the energy. 
In general then we choose to work with pseudospins in a (again uniformly) rotating reference 
frame with 




and m z . We obtain 
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dx 



(6.19) 



To be consistent with the linearization assumptions, P must be chosen so that ^ E integrates 
to and this term can be dropped. 
In Fourier space 
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m y (-p + (P- Q)x) + rh y (-p + (Q - P)x) 
m 2 (-p + (P - Q)x) + m z (-p + (Q- P)x 
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y/AJ2{m y ((P ~ Q)x) - m y ((Q - P)f)}. 



(6.20) 



We first need to minimize this functional to determine the ground state. Since there is 
no term linear in m z it follows that m z (q) = in the ground state. For ground state 
calculations we could work with the XY model in both commensurate and incommensurate 
cases. Minimizing with respect to fh y we find that in the ground state m y depends only on 
x and that its x-dependence is determined in Fourier space by solving 



Qlm° y (q x ) + — (m° y (q x + Q s ) + m° y (q x - Q s 
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S, 



qx,—Qsl- 



(6.21) 



We restrict our attention here to results which are valid to leading order in t so that the 
second term on the left-hand-side of Eq. ( |6.21| ) can be dropped. This gives 



tn\[A 



2ip s Q 

which gives a ground state energy per unit area equal to 



(6.22) 
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(6.23) 



A A Ps (Q-Py ' 2 
Minimizing this energy with respect to P we find that the overall ground state occurs for 



P = 1 (tn) 2 = t^(Qc 
Q 2 p s Q 4 512 I Q 



(6.24) 
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This result agrees with the analytic expression given in Ref.[ [34]] after making several cor- 
rections of some typos in Eq. (2.7) and (2.8). For large parallel fields P approaches 0, Q s 
approaches Q and the ground state of the system asymptotically approaches the ground 
state in the absence of interlayer tunneling. We emphasize that as Q approaches Q c from 
above rh° can become large and the linearization approximations will fail even in the ground 
state. 

Having identified the ground state pseudospin functional we are able to calculate the 
collective mode energies and response functions. To leading order in t we find that 

Att£ 2 dE 



m y (q) 



h dm z (— q 
2 



h z (q) + (2(3 + p s q 2 - p s P 2 ) m z (q) + y [m z (q + Q s x) + m 2 (q - Q s x )] 
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p s q 2 5m y (q) + — [Srh y (q + Q s x) + 8rh y (q - Q s x)} 
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The soliton lattice acts like an internal gratingEil which couples collective excitations of the 
ground state in the absence of interlayer tunneling whose wavevectors are separated by 
multiples of Q s ~ Q. The collective modes of the system for small t will consist of the zone- 
folded modes of the t = system with small corrections due to mode-mode coupling. The 
response functions of the system are readily evaluated numerically from the above equations. 
Artificial external gratings are often@ used to allow the infrared light to couple to finite- 
wavevector excitations of two-dimensional electron systems. The soliton lattice appears to 
offer an opportunity to couple to the Goldstone collective mode of the t — system at 
a wavevector which can be tuned by the application of an in-plane magnetic field. We 
reemphasize that as Q approaches Q c from above rh° can become large and the linearization 
approximations will fail even in the ground state. 
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VII. EXACT DIAGONALIZATION STUDIES 

We now turn to a discussion of exact numerical studies of the commensurate- 
incommensurate transition, and in particular the Bu dependence of the energy gap. The 
importance of long length scales in this transition will limit the power of the exact diagonal- 
ization approach, but we are still able to obtain some useful insights into the experimental 
results obtained by Murphy et al.@. We simulate a bilayer 2D quantum Hall system in the 
presence of a tilted magnetic field with a finite number of electrons, and perform our calcu- 
lations on the torus. We make use of the formalism developed by Haldane^l for the torus 
to block-diagonalize the Hamiltonian in each momentum sector. As noted in Sec. 0, the in- 
plane component of the magnetic field enters only in the tunneling matrix elements. In order 
to keep the momentum a good quantum number, we are constrained to use only discrete 
values of the magnetic field, namely the values which correspond to an integer number of 
flux quanta enclosed between the layers. We therefore study the transition for a fixed value 
of B» , and continuously vary the tunneling between the layers. This is to be contrasted with 
the experimental situation in which the tunneling is fixed and Bn varies continuously. 
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Fig. (|) shows the pseudo-magnetization calculated as a function of tunneling amplitude 
for various fixed values of B», for ten electrons and d/£ = 1. For fixed B\\ the commensurate 
state which optimizes the tunneling energy at the cost of exchange energy occurs at large 
t. We see in Fig. (§) that the component of pseudo-magnetization aligned with the effective 
Zeeman field increases with increasing t and decreases with B\\, even in the commensurate 
state. This quantum fluctuation effect is not captured in the classical field theoretic results 
we have presented since, , strictly speaking, they apply only when the tunneling amplitude is 
small. The increase in the effective magnitude of the ordered moment with tunneling in the 
commensurate state may explain the discrepancy between the t 1 ^ 2 behavior of the ordered 
moment predicted here and the approximately linear behavior seen experimentally.B 

In our finite-size studies the phase transition to the incommensurate state appears as a 
level crossing between states with different Haldane pseudomomenta which is accompanied 
by a large decrease in the component of the pseudospin aligned with the Zeeman field. In 
the thermodynamic limit and close to the phase transition, the pseudospin in the incom- 
mensurate state is expected to be aligned with the Zeeman field except in the domain walls 
of the soliton lattice. Thus the spatial average of aligned moment is expected to decrease 
continuously at the phase boundary. 

In our finite size calculations the wavevector Q corresponding to the rotating Zeeman 
field satisfies, 

/27T\V2 

^ = ^(-j (7.1) 

where N is the number of states in a Landau level, = B\\dL/$> , and L is the edge length 
for a square finite size square sample. For N = 10 Q£ ~ 0.79 even for = 1. Even 
the smallest value of Q£ that we can consider is close to the edge of the regime where we 
expect the long wavelength theory developed in previous sections to apply. We do however 
still expect the phase transition to occur outside this regime. In order to compare these 
microscopic calculations more directly with the long wavelength theory we have evaluated 
the energy gap as a function of t. We present our results in terms of an estimate of the 
dimensionless parameter Q£ which plays the central role in the long- wavelength theory. 

QZ = — (—r 1 ) ■ ( 7 - 2 ) 



To plot our results we have evaluated in the Hartree-Fock approximation; t = Asas/2 
corresponding to full pseudospin polarization and 

1 POO 

where = (27re 2 /e/c) exp(— kd), and h(k) = — exp(— k 2 /2). Both p s and t are substantially 
reduced by quantum fluctuations. Previous estimates from finite-size exact diagonalization 
calculations^ suggest that t is reduced by a larger fraction than p s , so that the Hartree- 
Fock approximation should increasingly underestimate Q£ as the layer separation increases. 
According to the P-T model, the transition should occur when Q£ = 4/tt. This will allow 
us to compare the HF prediction with the exact diagonalization result. 
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Figs. H B an d § show the Q^hf dependence of the energy gap for d/£ = 0.5, 1 and 2 
respectively. Note also that the system studied by Murphy et al. corresponds approximately 
to d/£ = 1.85. The cusp in the energy gap seen in these figures results from the ground state 
level crossing and is not obviously associated with the square root singularities in the energy 
gap predicted by Read0 which would be expected to appear only in much larger systems 
than we are able to treat. Nevertheless it does show the expected behavior of weakening 
at larger layer separations where the spin-stiffness and hence the meron-pair string energy 
is weakened. We also notice that Hartree-Fock theory accurately predicts the critical field 
for d/£ — 0.5, becomes worse at djl = 1 and fails badly at d/£ = 2. The direction of 
the discrepancy is in the direction anticipated by the above discussion since the value of 
Q^hf at the transition is smaller than 4/tt. These results suggest that the mean source of 
discrepancy between the experiments of Murphy et al. and the P-T theory with Hartree- 
Fock parameters is quantum fluctuations which reduce the ordered moment aligned by the 
effective magnetic field. This reduction is responsible for a proportional reduction in the 
energy gained by forming the commensurate state and reduces the parallel field strength at 
which the transition occurs. For parameters appropriate to the system studied by Murphy 
et al. our P-T theory with p s and t given by Hartree-Fock theory gives B» = 1.3T whereas 
the experimental value is BT, « 0.8T. From Fig. ([S]), we find B*\ w 0.6 T, in substantially 
better agreement with experiment. 

VIII. SUMMARY 

We have presented in this paper a theory of the interesting effects of a weak interlayer 
tunneling in double layer quantum Hall systems that spontaneously develop interlayer co- 
herence in the absence of tunneling. We have discussed the properties of the ground state, as 
well as low energy collective excitations (neutral and charged) of the system, using both effec- 
tive field theory and microscopic Hartree-Fock approaches. In particular, we have identified 
a novel commensurate-incommensurate phase transition driven by an in-plane component 
of the external magnetic field, which has been observed in recent experiments. Our theory 
is in good qualitative and semiquantitative agreement with experiments. 
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TABLES 



TABLE I. Charged Spin Texture Energies at vt = 1 for Double Layers Systems with Tunneling. 
p s = p s j {e 2 /£) and t = t/(e 2 /£) where p s is the pseudospin stiffness, t is the renormalized tunneling 

amplitude, £ is the magnetic length, To = 8p s /£ is the soliton string tension and £ = ^ 27r ^ p« ^ ^ i s 
the domain wall width. 
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FIGURES 



FIG. 1. Schematic conduction band edge profile for a double-layer two-dimensional electron 
gas system. 

FIG. 2. Phase diagram for the double layer QHE system (after Murphy et al.). Only samples 
whose parameters lie below the dashed line exhibit a quantized Hall plateau and excitation gap. 

FIG. 3. A process in double-layer two-dimensional electron gas systems which encloses flux 
from the parallel component of the magnetic field. The quantum amplitude for such paths is 
sensitive to the parallel component of the field. 

FIG. 4. Illustration of a meron pair in the presence of tunneling which confines the region of 
spin twist to a relatively narrow domain wall or 'string'. Each end of the string is a vortex carrying 
charge ±1/2. 

FIG. 5. Magnetization as a function of tunneling, for several integer numbers of flux quanta 
between the layers at iV = 10 and djt = 1. The curves are labeled by the integer number of flux 
quanta produced by the parallel component of the field in the finite size system. 

FIG. 6. The energy gap as a function of Q£hf, for By corresponding to one enclosed flux 
quantum, i.e. = 1, for N = 8 and djl = 0.5. 

FIG. 7. The energy gap as a function of Q£hf, for B\\ corresponding to one enclosed flux 
quantum, i.e = 1, for N = 8 and d/l = 1. 

FIG. 8. The energy gap as a function of Q£,hf, for B\\ corresponding to one enclosed flux 
quantum, i.e = 1, for TV = 8 and d/i = 2. 
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